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摘 X. 黄河 流域 地 域 广阔 ,但 气象 站 点 分 布 较 少 ,导致 气象 资料 短缺 。 卫 星 降 水 可 以 作为 气象 站 
点 观测 的 重要 补充 ,但 其 空间 分 辨 率 有 限 ,导致 其 在 区 域 研究 中 作用 有 限 。 以 黄河 流域 作为 研究 


区 域 , 针 对 全 球 降 水 观测 计 靶 
年 、 标 准 年 以 及 温润 年 3 个 典型 气候 年 份 ,在 


(GPM) 卫 星 降 水 产品 ,以 2002、2012 年 和 2020 年 降水 数据 作为 干旱 


综合 考虑 归 一 化 植被 指数 (NDVI1) 、 数 字 高 程 模 型 


ana 


(DEM) ,坡度 (Slope)、 地 表 温 度 (LST) 和 风速 (WDS) 多 种 反映 降水 量 空间 分 布 特征 的 影响 因子 及 其 
空间 非 平 稳 性 特征 的 基础 上 ,采用 地 理 加 权 回 归 (GWR ) 模 型 .混合 地 理 加 权 回 归 (MGWR ) 模 型 2 


种 降 尺 度 方法 ,得 到 了 黄河 流域 1km 空间 分 辨 率 的 降 尺 度 降水 数据 ,并 进一步 通过 地 鱼 
数据 对 降 尺 度 结果 进行 验证 。 结 果 表 明 :(1) GPM 年 降水 数据 与 地 本 
2012 年 和 2020 年 的 黄河 流域 地 区 具有 较 高 的 相关 性 。(2) 
辨 率 得 到 了 显著 提高 , 且 在 降水 变化 的 空间 细节 表达 方 玫 


气象 站 点 
气象 站 点 观测 数据 在 2002、 
经 MGWR 模型 降 上 尺度 的 降水 数据 空间 分 
较 GWR 模型 更 优 ,(3) 在 3 个 典型 气候 


年 份 中 ,MGWR 模型 在 降水 量 标准 年 中 相对 于 GWR 模 型 具有 更 高 的 准确 性 。 研 究 结果 能 够 为 相 
关 区 域 范 围 的 降水 降 尺 度 研究 提供 宏观 参考 与 借鉴 ,促进 区 域 气候 水 文 研究 。 
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降水 作为 全 球 水 循环 过 程 中 的 关键 要 素 之 一 ， 
获取 高 时 空 分 辨 率 的 降水 数据 在 生态 学 ,水 文学 和 
气象 学 研究 等 方面 具有 重要 意义 。 传 统 观 测 降 
水 方法 主要 包括 雨量 计 和 雷达 2 种 ,然而 这 类 观测 
方法 站 点 数量 有 限 且 空间 分 布 不 均 , 使 得 高 空间 分 
辨 率 的 降水 数据 资料 匡 乏 ,在 研究 降水 方面 数据 受 
到 限制 ,同时 ,降水 具有 高 时 空 异 质 性 ,传统 观 
测 降 水 方法 难以 准确 反映 时 空 分 布 规律 。 随 着 卫 
星 技术 的 发 展 ,遥感 反 演 降水 为 获取 降水 资料 提供 
了 一 种 新 途径 。 卫 星 和 遥感 数据 具有 不 受 天 气 和 复 
杂 地 形 影 响 、 陆 地 海洋 全 履 盖 大 面积 全 天 候 观测 
等 特点 ,因而 成 为 获取 降水 时 空 变化 信息 的 重要 方 
AO, 
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近 几 十 年 来 ,许多 机 构 和 政府 组 织 开 发 了 一 
系列 区 域 和 全 球 的 降水 数据 集 , 例 如 :利用 人 工 神 
经 网 络 从 遥感 卫星 产品 中 估算 降水 量 ;全球 降 
水 气候 学 项 目 (CGPCP) ;热带 降水 观测 任务 
(TRMM) 与 全 球 降水 观测 计划 (CPMJ)5 等 产品 ， 
这 些 降 水 产品 为 相关 研究 提供 了 有 力 的 数据 文 
TR. GPM 卫星 计划 是 由 美国 宇航 局 (NASA ) 执 行 ， 
日 本 宇宙 航空 开发 机 构 (JAXA) 以 及 欧盟 等 其 他 一 
些 国际 组 织 合作 建设 的 主要 用 于 气象 预报 的 卫星 
监测 系统 。GPM IMERG (Integrated multi-satellit E 
retrievals for GPM) 降 水 产品 自发 布 以 来 ,已 经 在 世 
界 各 地 开展 了 大 量 GPM 卫星 产品 的 适用 性 研究 。 
Anjum 等 在 巴基斯坦 北部 进行 了 4 个 季节 的 
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TRMM 5 GPM 降水 数据 的 精度 评估 ,得 到 CPM 降水 
数据 的 降水 量 估算 更 加 可 靠 的 结论 ; Tan 等 在 新 
加 坡地 区 进行 了 TMPA (TRMM multisatellite precipi- 
tation analysis) Ej GPM 降水 数据 对 比 研 究 , 证 明 
GPM 降水 数据 在 降水 空间 变异 性 表现 与 降水 估算 
能 力 方面 具有 更 好 的 性 能 ;Chen 等 “的 研究 中 ， 
GPM 降水 数据 可 以 更 加 精确 地 估算 中 国 降水 ,并 
且 在 新 疆 与 青藏 高 原 地 区 的 降水 量 估算 精度 显著 
提高 。 根 据 现 有 研究 内 容 ,GPM 降水 数据 在 不 同 
地 区 时 间 尺 度 和 地 形 条件 下 的 性 能 差异 很 大 。 
尽管 GPM 与 TRMM 降水 数据 相 比 ,提高 了 对 弱 降 水 
(«0.5 mm'h 和 固态 降水 的 探测 能 力 ,及 其 时 空 分 
PERE .观测 精度 和 准确 性 “ ,提供 了 更 加 可 靠 的 降 
水 估计 ,但 是 它们 的 空间 分 辩 率 往往 过 于 粗糙 ,无 
法 呈现 降水 空间 分 布 的 细微 变化 ,所 以 有 必要 对 卫 
星 降 水 数据 产品 进行 降 尺 度 处 理 , 以 提高 其 精度 和 
空间 分 辩 率 ” ,服务 于 区 域 尺度 或 更 小 尺度 的 精 
细 化 水 文 分 析 。 

降 尺 度 方法 大 体 上 可 以 分 为 动力 降 尺 度 和 统 
计 降 尺度 两 大 类 ,相对 于 动力 降 尺 度 方法 ,统计 
降 尺 度 因 其 具有 计算 量 小 节省 时 间 方法 众多 形式 
灵活 .更 易于 操作 等 优点 而 被 广泛 应 用 于 卫星 遥感 
降水 产品 的 降 尺 度 人 研究 中 。 地 理 加 权 回 归 (GWR ) 
模型 是 具有 代表 性 的 统计 降 尺 度 方法 ,被 广泛 应 用 
在 卫星 降水 数据 产品 的 降 尺 度 研究 中 。Immerzeel 
等 "通过 建立 TRMM 降水 数据 与 NDVI 之 间 的 指 
数 模型 ,获取 了 伊比 利 亚 半岛 范围 内 空间 分 辨 率 为 
1 km 的 年 降水 数据 ;Jia 等 "在 此 基础 上 ,通过 构建 
TRMM 数据 与 归 一 化 植被 指数 (NDVI) 、 数 字 高 程 
模型 (DEM ) 的 多 元 线性 回归 模型 使 柴 达 木 贫 地 内 
降水 数据 的 空间 分 辨 率 提高 至 1 km; 温 伯 清 等 
SEF GWR 模型 对 陕西 地 区 的 GPM 降水 数据 进行 
降 尺 度 处 理 , 发 现 GWR 降 尺 度 降 水 数据 在 陕西 省 
内 精度 较 好 ,能 够 较为 准确 地 反映 陕西 省 内 的 降水 
分 布 ; 曾 昭 昭 等 3 采用 GWR 模型 进行 秦 巴 山区 
TRMM 降水 数据 降 尺 度 研 究 ,发 现 GWR 降 尺 度 降 
水 数据 具有 较 强 的 适用 性 ;内 路 明 等 “基于 CWR 
模型 对 中 国 三 大 流域 的 TRMM 和 GPM 降水 数据 进 
行 降 尺度 处 理 , 得 到 了 三 大 流域 1 km 分 辨 率 的 降 
尺度 降水 数据 。GWR 模型 适合 在 环境 变量 具有 明 
显 的 空间 非 平 稳 性 下 进行 建 模 ,但 是 在 许多 实际 问 
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题 中 ,环境 变量 包括 全 局 变量 和 局 部 变量 2 种 类 型 , 
并 且 环 境 变 量 存在 不 具有 空间 非 平稳 性 或 者 空间 
非 平稳 性 较 小 可 忽略 不 计 的 情况 。 降 水 是 一 种 复 
杂 的 自然 现象 ,其 特点 是 在 时 间 和 空间 上 都 具有 显 
著 的 可 变性 ,并且 影 响 降水 的 因素 中 可 能 同时 
包含 空间 平稳 和 空间 非 平 稳 2 种 类 型 。Arshad 577 
在 印度 河 倪 地 通过 建立 TRMM 降水 数据 与 NDVI、 
DEM ,坡度 (Slope) .地表 温度 (LST) 和 风速 (WDS ) 
AY GWR 与 混合 地 理 加 权 回 归 (MGWR ) 模 型 ,并 通 
过 地 理 比 率 分 析 (GRA) 和 地 理 差 异 分 析 (GDA ) 将 
降 尺 度 降 水 数据 与 雨量 站 点 降水 数据 进行 合并 , 进 
一 步 提高 了 降 尺 度 降 水 数据 的 准确 性 ,因此 ,本 文 
引入 MGWR 模型 对 卫星 降水 产品 进行 降 尺 度 效果 
研究 。 

黄河 流域 是 我 国 的 三 大 流域 之 一 也 是 降水 的 
主要 承载 地 。 因 此 ,本文 以 黄河 流域 作为 研究 区 
域 , 以 地 面 气象 观测 站 点 降水 数据 作为 真实 值 , 通 
过 与 主流 的 GWR 模型 降 尺 度 结果 以 及 原始 的 GPM 
降水 数据 进行 对 比 ,分 析 MGWR 模型 在 不 同 干 湿 气 
候 条 件 下 对 黄河 流域 CPM 降水 数据 的 降 尺度 效果 ， 
以 期 为 相关 区 域 范围 的 降水 降 尺 度 研 究 提 供 宏观 
参考 与 借鉴 ,促进 区 域 气候 水 文 研究 。 


1 研究 区 概况 


本 文 以 黄河 流域 (95°50’~119°06'E,32°06'~ 
41?48'N ) 作 为 研究 区 域 , 黄 河 发 源 于 巴 彦 喀 拉 山脉 ， 
依次 流 经 青藏 高 原 .内 蒙古 高 原 和 黄土 高 原 、 下游 
冲积 平原 3 级 阶梯 (图 1)。 黄 河流 域 属于 我 国 典型 
的 干旱 、 半 干旱 、 半 湿润 气候 区 ,流域 总 面积 约 
79.46x10* km 。 根 据 黄河 水 利 委 员 会 提供 的 划分 方 
案 将 黄河 流域 划分 为 上 中 下游 3 个 区 域 ”。 

黄河 流域 地 域 广阔 ,整体 地 势 西高 东 低 ,起 伏 
显著 且 地 貌 复杂 ,由 西向 东 依 次 分 为 西部 、 中 部 、 东 
部 ,3 个 地 区 分 别 以 青藏 高 原 区 .黄土 高 原 区 和 黄 淮 
海平 原 区 为 主 。 由 于 黄河 流域 受 大 气 环流 和 季风 
环流 影响 的 情况 比较 复杂 ,流域 内 东西 气候 差异 显 
著 ,降水 空间 分 布 不 均 光 ,但 是 流域 蒸发 能 力 强 , 造 
成 流域 内 旱灾 事件 影响 面积 广 .持续 时 间 长 并 且 频 
繁 发 生 ”™ 1。 因此 ,获得 精确 、 快 速 的 降水 数据 对 黄 
河流 域 的 干旱 监测 至 关 重 要 。 
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图 例 
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注 :该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 的 审 图 号 为 CSC022)4312 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 
图 1 黄河 流域 气象 站 点 分 布 


Fig. 1 Distribution of meteorological stations in the Yellow River Basin 


2 数据 与 方法 


2.1 数据 来 源 与 处 理 
2.1.1 降水 数据 本 文采 用 的 GPM 卫星 降水 数据 
(简称 GPM 数据 ) 来 源 于 NASA 网 站 (https://www. 
nasa.gov/)。GPM 卫星 计划 根据 其 数据 反 演 算法 不 
同 提供 了 4 个 级 别 的 产品 ,本 文选 取 IMERG 产品 的 
最 新 版 本 IMERG V06 (IMERG. final) P 2001—2020 
年 逐日 的 空间 分 辨 率 为 0.1°x0.1° 的 GCPM 数据。 该 
数据 于 2019 年 3 月 发 布 , 它 整合 了 可 用 的 GPM, 
TMPA 遥感 大 获取 的 信息 ,并 进行 互相 校准 、 合 并 、 
插值 和 融合 生成 一 致 的 降水 数据 ,实现 了 2000 年 6 
月 至 今 的 降水 估 测 后 。 

根据 CPM 数据 提取 黄河 流域 2001 一 2020 年 的 
年 累积 降水 量 信息 ,计算 每 年 的 年 平均 降水 量 及 
2001 一 2020 年 平均 降水 量 ,通过 比较 年 平均 降水 量 
5j 2001—2020 年 平均 降水 量 , 选 取 3 个 典型 气候 年 
份 (以 降水 量 值 最 低 的 2002 年 作为 干旱 年 ,与 年 平 
均 降 水 量 值 最 接近 的 2012 年 作为 标准 年 ,降水 量 值 
最 高 的 2020 年 作为 湿润 年 ) ,进行 黄河 流域 GPM 数 
据 在 不 同 降 尺 度 模型 下 的 降 尺 度 精 度 研 究 (图 2)。 
2.1.2 环境 变量 本文 NDVI 与 LST 均 来 源 于 NASA 
(https://ladsweb.modaps.eosdis.nasa.gov/) , 其 中 NDVI 
采用 的 是 MOD13A3 月 合成 植被 指数 数据 ,空间 分 辨 
率 为 1 km; MODIS LST 系列 数据 产品 中 MOD11A2 


LST, 时 间 分 辨 率 为 8 d, 空 间 分 辨 率 为 1 km。 运 用 
MRT (Modis reprojection tool) 软件 ,对 MOD13A3 与 
MOD11A2 数 据 进行 预 处 理 ,NDVI 采 用 最 大 值 合 成 
法 生成 。SRTM DEM 数据 来 源 于 地 理 空间 数据 云 
(https://www.gscloud.cn/) ,空间 分 辩 率 为 90 m, 运 用 
ArcGIS 进一步 获取 坡度 (Slope) 数 据 。 

2.1.3 站 点 数据 气象 站 点 观测 数据 来 源 于 国家 气 
象 信 息 中 心 (http://data.cma.cn) ,本 文选 取 黄 河流 域 
内 86 个 气象 站 点 的 降水 数据 以 及 WDS 数 据 ,其 中 
通过 对 WDS 数 据 进行 插值 处 理 来 获取 双 DS 栅 格 数 
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图 2 2001 一 2020 年 黄河 流域 年 降水 累积 量 
Fig. 2 Annual precipitation accumulation in the Yellow River 
Basin from 2001 to 2020 
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据 作为 环境 变量 ,降水 数据 作为 真实 值 对 降 尺度 前 
后 的 GPM 卫星 降水 产品 进行 精度 验证 。 
2.2 研究 方法 
2.2.1 降 尺 度 方法 

(1) GWR 模型 

GWR 模型 是 Brunsdon 等 于 1996 年 提出 的 用 
来 探测 空间 非 平稳 性 的 一 种 局 部 回归 分 析 方 法 ,是 


与 局 部 变量 。 本 文 分 别 对 2002、2012 年 和 2020 年 
进行 地 理 谈 异性 测试 检验 , 当 差 异 标准 值 小 于 0 时 ， 


说 明 辅 助 变量 的 空间 变异 性 较 弱 , 则 将 该 变量 作为 
全 局 变量 ;反之 , 则 将 该 变量 作为 局 部 变量 ,具体 如 


表 1 所 示 。 
2.2.2 降 尺 度 流 程 ” 降 尺度 方法 基于 2 个 基本 假设 
进行 :首先 ,降水 量 与 环境 变量 之 间 存 在 可 以 用 回 


传统 线性 回归 模型 的 扩展 。 其 基本 思想 是 变量 间 
的 关系 随 着 空间 位 置 的 变化 而 变化 ,在 回归 中 通过 
引入 空间 关系 权重 ,将 全 局 模型 中 的 参数 设 定 为 地 
理 位 置 的 函数 ,用 以 估算 研究 区 内 任意 位 置 的 相关 
变量 与 解释 变量 的 参数 在 空间 中 的 变化 。 已 有 研 
究 表明 ,NDVI.DEM LST Slope\WDS 是 影响 降水 的 
重要 因素 ,因此 本 研究 利用 这 5 个 环境 变量 作为 辅 
助 因 子 与 GPM 建立 GWR 模型 ( 式 1)。 
Y,=f,u,,0,) + Y Ao), t e(u,,v,) 
i-1,2,:,n 
式 中 :了 为 第 i 点 的 降水 量 (mm) ; (us) 29 98 i DEL 
点 的 经 纬度 坐标 ;p(wis0) 为 第 数 项 ;p(w 为 第 i 个 位 
置 点 上 的 第 i 个 回归 参数 ;大 ;为 第 k 个 空间 因素 在 第 
i 个 位 置 点 上 的 值 ;e(wisy) 为 每 个 位 置 上 的 残 差 。 
(2) MGWR 模型 
1999 年 Brunsdon 等 在 GWR 模型 基础 上 提 
出 了 MGWR 模 型。 在 实际 问题 中 ,空间 关系 中 存 
在 对 因 变 量 影 响 具 有 全 局 性 质 的 解释 变量 ,同时 存 
在 影响 是 局 部 的 解释 变量 ,因此 对 一 个 地 理 问 题 进 
行 完 整 的 空间 建 模 需 要 在 模型 中 进行 全 局 变 


(1) 


局 部 变量 的 区 分 。 降 水 具有 明显 的 空间 异 质 性 ， 
响 降水 的 因素 可 能 同时 存在 空间 平稳 与 空 "—— 


稳 2 种 类 型 ,因此 MGWR 模型 通过 引入 全 局 变量 ， 
是 能 够 解决 空间 非 平稳 性 问题 的 有 效 数据 分 析 方 
法 ( 式 2)。 


Y= 全 局 变 zx x. 


局 部 2 T» P(u, v )X, + elu; ») 


SUP X Xi BAB j LA RI SER YEAR i AEA 
EE ; œ 2 4 Jeg FY E83 RR S Bis) A Jr GL] 
归 的 回归 参数 ;j=1, 2, ++, k,l=1,2, ++, po 

本 人 研究 中 ,通过 地 理 变 异性 测试 计算 差异 标准 
值 ,根据 差异 标准 值 的 正 负 区 分 研究 区 内 全 局 变量 


(2) 


归 模 型 模拟 与 处 理 的 空间 关系 ;其 次 ,在 低空 间 分 
辨 率 下 建立 的 模型 可 以 利用 环境 变量 数据 进行 高 


分 辩 率 下 的 降水 预测 。 本 文通 过 MCWR 模型 与 

GWR 模型 对 GPM 数据 进行 降 尺 度 操作 ,图 3、 图 4 

中 的 流程 图 总 结 了 MGWR 和 GWR 模型 降 尺 度 流 
程 ,具体 操作 步骤 如 下 : 


表 1 辅助 因子 的 差异 标准 值 


Tab.1 Difference standard values of accessory factors 


辅助 因子 2002 年 2012 年 2020 年 
DEM -13792.4780 111.5426 -38.1083 
NDVI -10307.8175 -14268.9533 -11486.5637 
WDS -13864.9890 -12589.8497 -11786.7153 
Slope 283.4418 329.9034 374.2292 
LST 291.0240 -6264.2374 227.1289 


IE: DEM 为 数字 高 程 模型 ;NDVI 为 归 一 化 植被 指数 ;WDS 为 风速 ; 
Slope 为 坡度 ;LST 为 地 表 温 度 。 下 同 。 


GPM 
(0.19x0.19) 


地 理 变 异性 测试 


(0.19x0.19) (1 kmx1 km) 


全 局 变量 与 局 部 变量 
1 km 预测 降水 数据 


注 :DEM 为 数字 高 程 模型 ;NDVI 为 归 一 化 植被 指数 ;WDS 为 风速 ;Slope 
为 坡度 ;LST 为 地 表 温 度 ;GPM 为 全 球 降水 观测 计划 ;MGWR 为 混合 
地 理 加 权 回 归 。 下 同 。 

图 3 MGWR 模 型 降 尺 度 流 程 图 
Fig. 3 Downscaling flow chart of MGWR model 


= 


202308.00020v1 


chinaXiv 


ChinaX iv EAT 


1056 T 4 w 4635 


1 km 预测 降水 数据 


注 :GWR 为 地 理 加 权 回 归 。 下 同 。 
图 4 GWR 模 型 降 尺 度 流程 图 
Fig. 4 Downscaling flow chart of GWR model 


(1) ¥ NDVI,LST, DEM ,Slope .WDS 数 据 采 用 
克 里 金 插值 方法 重 采 样 至 与 CPM 数据 相同 的 空间 
分 辨 率 , 以 进行 匹配 。 

(2) 通过 地 理 变 异性 测试 筛选 出 环境 变量 中 的 
全 局 变量 与 局 部 变量 ,并 建立 重 采样 后 环境 变量 与 
GPM 数 据 之 间 的 关系 。 

(3) 将 重 采 样 后 环境 变量 数据 与 原始 GPM 数 
据 建立 GWR 模型 ,得 到 GWR 模 型 常数 项 .环境 变量 
对 应 系数 项 以 及 残 差 结果 ;根据 全 局 变量 与 局 部 变 
量 的 划分 ,将 重 采 样 后 全 局 变量 与 局 部 变量 按照 式 
(2) 建 立 MGWR 模型 ,得 到 MGWR 模型 的 常数 项 、 
环境 变量 对 应 系数 项 以 及 残 差 结 

(4) 将 0.1?x0.1 空 间 分 辩 率 的 CWR MGWR 模 
型 的 常数 项 .系数 项 以 及 残 差 数据 进行 栅 格 化 处 
理 , 并 运用 克 里 金 插值 法 重 采样 至 1 km. 

(5) 根据 式 (1) 与 式 (2) 分 别 进行 计算 ,得 到 
GWR 5 MGWR 模型 降 尺 度 处 理 后 空间 分 辩 率 为 
1 km 的 预测 降水 数据 。 

2.2.3 评价 指标 为 了 检验 降 尺 度 结 果 的 精度 ,将 
地 面 气象 站 点 实测 数据 作为 真实 值 与 降 尺 度 结果 
进行 比较 验证 。 采 用 决定 系数 ( 尼 ) 偏差 (Bias ) 以 
及 均 方 根 误差 (RMSE)3 个 指标 作为 评价 指标 。 尼 又 
PUA DERE ,用 以 揭示 降 尺 度 结果 与 真实 值 之 间 的 
相关 程度 , 值 越 大 表明 两 者 之 间 相 关 性 越 高 。Bias 
表示 降 尺 度 结果 与 真实 值 之 间 的 偏离 程度 , 值 越 


小 ,表明 两 者 之 间 差 异 越 小 。RMSE 用 以 表示 降下 
度 结果 和 降水 真实 值 的 偶 差 ,其 值 越 小 代表 两 者 越 


接近 。 

VE- DY,-D) 

R = = (3 ) 
[ 2.0 - Xy, - YY 
Bias = = -1 (4) 
ix 

y= 

RMSE= |=! " (5) 


AP Xu X DS ss BOSE DM S ES BE Je: GPM 降水 数 
A MH EC OE TR LY. Y 为 气象 站 点 数据 FE; n 
为 气象 站 点 数 。 


3 结果 与 分 析 


3.1 原始 数据 精度 验证 

GPM 数据 在 黄河 流域 内 降水 空间 分 布 如 图 
Sa-c 所 示 。 利 用 研究 区 内 86 个 地 面 气 象 观测 站 点 
降水 量 值 为 真 值 对 其 进行 验证 (图 5d~f)。 在 黄河 
流域 2002 .2012 年 和 2022 年 3 个 典型 气候 年 中 CPM 
数据 的 尼 均 大 于 0.7; Bias 4} 3l] X 0.026, 0.300, 
0.311; RMSE 433] 23 62.144 .72.644 .115.091。 总 体 
上 ,CPM 数据 与 真实 值 存在 较 好 的 相关 性 ,但 是 由 
于 降水 的 非 线 性 特征 、 星 载 雷 达 的 探测 性 能 和 反 
演技 术 的 局 限 性 ,CPM 数据 和 实际 降水 存在 一 定 
差异 。 
32 降 尺 度 结果 及 其 精度 分 析 
3.2.1 GWR 模 型 降 尺 度 结果 及 其 精度 分 析 采用 
GWR 模型 对 GPM 数据 进行 降 尺 度 处 理 , 得 到 了 其 
在 黄河 流域 2002.2012 年 和 2022 年 3 个 典型 气候 年 
份 中 的 1km 空 间 分 辩 率 预测 年 降水 数据 空间 分 布 
图 (图 6a~c)。 从 图 中 可 见 ,GWR 降 尺 度数 据 ( 以 下 
简称 GWR 数据 ) 空 间 分 状 率 有 明显 的 改善 ,从 原始 
的 0.1°x0.1° 提 升 至 1 kmx1 km; 同 时 ,降水 空间 分 布 
趋势 与 GPM 数据 整体 上 保持 一 致 ,并 且 随 着 空间 分 
辩 率 的 提高 ,降水 空间 分 布 细节 相 较 于 CPM 数据 表 
现 更 加 清晰 。 
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(b) 2012 年 GPM 原 始 数据 N 
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(c) 2020 年 GPM 原 始 数据 


图 例 降水 量 /mm 图 例 降水 量 /mm 图 例 降水 量 /mm 
0 200 km Cj 714.089 0 200 km ”器 黄河 边界 g 0 200km 号 黄河 边界 1519.080 
CO XA 169.962 口 省 区 边界 168.048 口 省 区 边界 168.795 
(d) 2002 年 GPM 原 始 数 据 精 度 验证 (e) 2012 年 GPM 原 始 数 据 精 度 验 证 (f) 2020 年 GPM 原 始 数据 精度 验证 
600r R2=0.700 900r p20 786 1600F R2-0.736 
550} Bias-0.006 * a. Me 800 - Bias=0.300 . E 1400} Bias=0.311 
E 500 - RMSE-62.144* -- E RMSE=72.644 = 1200 上 RMSE=115.091 
rpm | = VL. gg 700 - dü 
x EC * 600l X 1000F 
x 400 - B 
= 350L = 500+ = 800r 
Ay [as] A 
© 300} 9 soo} © 6p 
250+ 300 4007 = 
1 L 1 1 1 L 1 1 L J I L 1 f 


h 
RS 


气象 站 点 降水 量 /mm 


图 5 黄河 流域 GPM 卫星 年 降水 量 空 


aP i 1 i i 
OD HD DH OP PD DY L 
PSSST OCH P SHS 


气象 站 点 降水 量 /mm 


SS a SSS ee ae | 
N > o © eO 
SE EE SF SS 

气象 站 点 降水 量 /mm 


s 间 分 布 及 其 精度 验证 


Fig. 5 Spatial distribution and accuracy verification of GPM satellite annual precipitation in the Yellow River Basin 


(a) 2002 年 GWR 降 尺度 数据 À 


(b) 2012 年 GWR 降 尺度 数据 N 


A 


(c) 2020 年 GWR 降 尺度 数据 


图 例 降水 量 /mm 图 例 降水 量 /mm 图 例 降水 量 /mm 
0 2 200km 器 黄河 边界 gr 927.207 0 200km 黄河 边界 E 0 200km 器 黄 河 边界 | 
口 省 区 边界 M 3] 687 器 省 区 边界 154.756 中 省 区 边界 198.370 
(d) 2002 年 GWR 降 尺度 数据 精度 验证 (e) 2012 年 GWR 降 尺度 数据 精度 验证 (D2020 年 GWR 降 尺度 数据 精度 验证 
E 650r R20.670 : 900r 8-0.773 Jen R^-0.683 
600 Bias-0.029 & 800. Bias=0.300 = £1400} Bias=0.309 . 
hus RMSE-73.940 Ti RMSE-134.04 la 
Æ 700L ¥ 1200 
a 600 L z 1000 
Hi 500L Hi 300 
É 400. dM 
© O 
300. a 400r x 
eT i a jt Rp ee 
9 NJ IN] 
A S SIC C PRU SOS SS S ge n gv eg $ SS SS © gg SF SF sS 
气象 站 点 降水 量 /mm 气象 站 点 降水 量 /mm 气象 站 点 降水 量 /mm 


图 6 黄河 流域 GWR 降 


尺度 数据 年 降水 量 空 间 分 布 及 其 精度 验证 


Fig.6 Spatial distribution and accuracy verification of GWR downscaling data of annual precipitation in the Yellow River Basin 


通过 定量 分 析 GWR 数据 的 精度 (图 6d~f)。 可 年 ) 表 现 最 为 明显 ,降低 了 0.053; 而 在 标准 年 (2012 
以 看 出 ,在 3 个 典型 气候 年 中 ,GWR 数据 与 CPM 数 ” 年 ), 相 较 于 GPM 数据 及 仅 降 低 了 0.013, 同 时 在 
据 相 比 ,GWR 数据 相关 性 略微 降低 ,在 湿润 年 (2020 ”2012 年 中 CPM 数据 与 GCWR 数据 的 Bias 均 为 0.300， 
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WIH GWR 模型 在 降水 量 标准 年 表现 较 好 ; 降 尺 度 
后 3 个 典型 气候 年 份 的 RMSE 均 略微 提高 。 整 体 来 
看 ,GWR 模型 在 降水 量 标准 年 表现 效果 较 好 。 
3.2.2 MGWR 模型 降 尺 度 结 果 及 其 精度 分 析 K 
用 MGWR 模型 对 GPM 数据 降 尺 度 处 理 , 得 到 了 其 
在 黄河 流域 2002、2012 年 和 2022 年 3 个 典型 气候 年 
份 中 的 1 km 空间 分 辩 率 预测 年 降水 数据 空间 分 布 
图 (图 7a~c)。 从 图 中 可 见 ,MGWR 降 尺 度数 据 ( 以 
下 简称 MGWR 数据 ) 在 降水 空间 分 布 细 节 相 比较 于 
GPM 数 据 和 GWR 数据 有 明显 提升 ;同时 ,其 降水 空 
间 分 布 趋势 与 GPM 数据 整体 上 保持 一 致 ,并 且 与 
GWR 数据 相 比 ,MGWR 数据 降水 空间 分 布 细节 的 
表现 能 力 更 强 且 边缘 更 加 平滑 。 

通过 定量 分 析 MGWR 数据 的 精度 (图 7d~f)。 
可 以 看 出 ,在 3 个 典型 气候 年 中 ,MGWR 5 GWR Zt 
据 表 现 相同 的 是 ,在 黄河 流域 范围 内 降水 结果 的 相 
关 性 均 略 微 低 于 CPM 数据 ;表现 不 同 的 是 , MGWR 
数据 在 干旱 年 (2002 年 ) 尺 降低 表现 得 最 为 明显 , 比 
GPM 数据 降低 了 0.102; 而 在 湿润 年 (2020 年 ), MG- 
WR 数 据 尼 较 CWR 数据 有 略微 提升 ,由 0.683 提升 
至 0.688。MGWR 数据 在 3 个 典型 气候 年 份 的 RMSE 
均 略 微 有 提高 ,但 是 在 2020 年 2 种 降 尺 度数 据 的 


(a) 2002 年 MGWR 降 尺度 数据 


(b) 2012 年 MGWR 降 尺度 数据 


RMSE 值 差异 较 小 ,整体 上 来 看 ,说 明 MGCWR 模型 在 
降水 量 湿 润 年 表现 效果 较 好 。 
3.3 2 种 模型 降 尺 度 结果 对 比 

由 图 6a~c、 图 7a~c 可 知 ,GWR 5j MGWR 数据 ， 
空间 分 辩 率 均 从 0.1°x0.1° 提 升 至 1 kmx1 km; 并 且 
年 降水 空间 分 布 趋势 大 致 为 自 北 向 南 逐 渐 递 增 , 降 
水 量 峰 值 区 主要 位 于 黄河 流域 南端 ,在 标准 年 
(2012 年 ) 和 湿润 年 (2020 年 ) ,黄河 流域 东部 存在 另 
一 降水 量 峰值 区 ,这 与 CPM 数 据 空间 分 布 趋势 整体 
上 保持 一 致 ,并 且 随 着 空间 分 辩 率 的 提高 ,降水 空 
间 分 布 细节 相 较 于 CPM 数据 表现 得 更 加 清晰 ,一些 
小 范围 .局 部 降水 空间 分 布 细节 表现 得 更 加 明显 ， 
在 一 定 程度 上 ,可 以 更 加 清晰 地 反映 出 黄河 流域 降 
水 的 空间 分 布 规律 。 与 GWR 模型 相 比 ,MGWR 数 
据 空 间 细 节 表 现 能 力 更 强 ,降水 数据 空间 分 布 边缘 
更 加 平滑 ,同时 在 这 3 个 典型 气候 年 份 中 可 以 看 出 ， 
在 降水 量 较 丰 富 的 地 区 ,MGWR 数 据 能 够 更 好 地 表 
现 出 降水 量 空间 分 布 的 细节 效果 。3 个 典型 气候 年 
份 降 尺 度数 据 对 降水 空间 分 布 细节 表现 能 力 为 : 
MGWR 数据 >GWR 数据 >GPM 数据 ,MGWR 数据 可 
以 较 好 地 还 原 GPM 数据 在 黄河 流域 地 区 的 降水 观 
测 信息 ,并且 在 降水 量 丰 富 的 地 区 表现 更 好 。 


(c) 2020 年 MGWR 降 尺度 数据 


图 例 降水 量 /mm 图 例 降水 量 /mm 图 例 降水 量 /mm 
0 200km 加 黄河 边界 gar 773.758 0 200km TRIIR 894.014 0 200km 号 黄河 边界 gr 1378.530 
S 器 省 区 边界 ™ 144,850 器 省 区 边界 0.91917 — 口 省 区 边界 128.813 
(d) 2002 年 MGWR 降 尺度 数据 精度 验证 (e) 2012 年 MGWR 降 尺度 数据 精度 验证 (f) 2020 年 MGWR 降 尺度 数据 精度 验证 
E 600F 20,598 i 8 900 - 0.690 . £ 1400 - Ro0.688 
s; 5950F Bias-0.002  , " A -— 800 Bias- —0.092 gi 1200 Bias-0.154  。 
关 500} RMSE-71.950 = "ye" 34 700 + RMSE-92.269 X opp || RMSE=134.229 
8 450 "UP. " = 1000 s 
EN ot 600+ P 
T "Ud a a 500+ uM 
Š 350 : 400 - E 600 
E 300 & 300+ Z 400 : 
Š 250 Š M z Pu 
| 20 TD 200[L ea , ,; ; , | | j 200 4 1 | 1 1 f 
IN] S 
PPP PP FP PHP SS OS SFL EPCS SY SEF EE SF Ss 


气象 站 点 降水 量 /mm 


气象 站 点 降水 量 /mm 


气象 站 点 降水 量 /mm 


图 7 黄河 流域 MGWR 降 尺 度数 据 年 降水 量 空间 分 布 及 其 精度 验证 


Fig.7 Spatial distribution and accuracy verification of MGWR downscaling data of annual precipitation in the Yellow River Basin 
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通过 图 8 可 以 看 出 ,在 干旱 年 (2002 年 ) 与 标准 
年 (2012 年 ) 降 水 量 差 值 多 为 正 值 , 且 正 值 范围 分 布 
较 广 ;在 湿润 年 (2020 年 ) ,黄河 流域 南部 地 区 显示 
为 正 值 ,北部 地 区 负 值 相对 较 多 。 总 体 上 ,由 于 干 
旱 年 (2002 年 ) 总 降水 量 相对 更 小 ,其 差 值 的 值 域 范 
围 也 相对 更 小 ;湿润 年 (2020 年 ) 总 降水 量 相对 更 
大 ,其 差 值 的 值 域 范围 也 相对 更 大 。 原 因 可 能 是 : 
降水 量 的 多 少 有 影响 了 MCWR 5 GWR 数据 的 精度 ， 


(a) 2002 年 降 尺度 差 值 ， 


0 200 km 图 例 降水 量 差 值 /mm 0 200km 图 例 
CXXXI 一 - -369.742 


N (0) 2012 年 降 尺 度 差 值 。 


中 黄河 边界 
口 省 区 边界 


ChinaXiv 合 作 期 刊 


7 期 柏 和 荷 等 : 基于 MGWR 模 型 的 黄河 流域 GPM 卫星 降水 数据 降 尺 度 研究 1059 


而 导致 了 降 尺 度 结果 差 值 值 域 的 变化 。 

从 表 2 中 可 以 看 出 ,GWR MGWR 数 据 在 3 个 典 
型 气候 年 份 的 尼 值 均 介 于 0.5-0.8 之 间 ,标准 年 
(2012 年 ) 尼 最 高 ,干旱 年 (2002 年 ) 湿润 年 (2020 
年 ) 随 着 CPM 数据 尼 的 升 高 ,GCWR 5 MGWR 数据 的 
尼 也 随 之 有 一 定 的 提升 ,说明 MGWR 模型 在 降水 量 
湿润 年 表现 效果 较 好 ,其 降 尺 度数 据 精 度 对 原始 数 
据 的 精度 具有 一 定 的 依赖 性 。 


降水 量 差 值 /mm 0 200 km 图 例 降水 量 差 值 /mm 
- 291.443 —À ^ 口 黄河 边界 e 592.802 
" —600.420 Co IA —1193.730 


图 8 黄河 流域 降 尺 度 差 值 空间 分 布 


Fig. 8 Spatial distribution of downscaling differences in the Yellow River Basin 


3&2 不 同年 份 降 尺 度 结果 验证 
Tab. 2 Verification of downscaling results in 
different years 


评价 指标 年 份 GPM GWR MGWR 
R 2002 0.700 0.670 0.598 
2012 0.786 0.773 0.690 

2020 0.736 0.683 0.688 

Bias 2002 0.026 0.029 0.022 
2012 0.300 0.300 -0.092 

2020 0.311 0.309 0.154 

RMSE 2002 62.144 65.482 71.950 
2012 72.644 73.940 92.269 

2020 115.091 134.041 134.229 


注 : RAE AM Bias 为 偏差 ;RMSE 为 均 方 根 误差 ;GPM 为 全 
球 降水 观测 计划 ;MGWR 为 混合 地 理 加 权 回 归 ;GWR 为 地 理 加 权 
归 。 


在 3 个 典型 气候 年 份 中 ,GWR MGWR 数据 与 
GPM Ziff fH EE RC (535 FE fI , RMSE 升 高 ,说 明 GPM 
数据 与 真实 值 的 尼 并 没有 因 其 经 过 GWR 5j MGWR 
降 尺 度 处 理 而 提高 。 但 是 降 尺 度 结果 均 降低 了 BI- 
AS 值 ,说 明 降 尺度 结果 的 准确 性 得 到 了 提高 ,在 年 
尺度 上 利用 NDVI、LST、DEM Slope .WDS 对 GPM 数 
据 进 行 降 尺度 具有 一 定 的 科学 性 。 但 是 CWR 、MG- 
WR 数据 的 让 均 低 于 CPM 数据 ,其 原因 除了 实测 数 


据 在 一 定 程度 上 限制 了 降 尺 度 结果 的 验证 ,造成 一 
定 程 度 的 不 确定 性 之 外 ,还 可 能 是 由 于 黄河 流域 内 
降水 丰富 多 变 昌 地 势 复杂 导致 对 降水 的 反 演 较为 
困难 ,并 且 本 文中 引入 的 环境 变量 有 限 ,在 处 理 过 
程 中 存在 一 定 的 信息 丢失 ,导致 2 种 降 尺 度数 据 的 
尼 略 低 。 所 以 ,GWR .MGWR 数据 在 提升 分 辨 率 的 
同时 ,牺牲 掉 了 一 部 分 的 数据 精度 *。 


4 讨论 


目前 ,GWR 模 型 降 尺 度 方法 被 众多 学 者 广泛 应 
用 于 卫星 遥感 降水 产品 的 降 尺 度 研 究 中 3, 但 是 
其 使 用 受到 环境 变量 是 否 具 有 空间 平稳 性 的 限 
制 。 降 水 是 一 种 复杂 的 自然 现象 ,其 特点 是 在 时 间 
和 空间 上 都 具有 显著 的 可 变性 3, 并 且 影 响 降水 
的 因素 中 可 能 同时 包含 空间 平稳 和 空间 非 平稳 2 种 
类 型 。 因 此 ,分 析 MGWR 模型 在 不 同 干 湿 气 候 条 件 
下 对 黄河 流域 GCPM 数据 的 降 尺 度 效果 ,能 够 为 相关 
区 域 范围 的 降水 降 尺 度 研 究 提供 宏观 参考 与 借鉴 ， 
促进 区 域 气候 水 文 研究 。 

整体 上 ,本 文 证 实 了 MGWR 模型 相对 于 传统 
GWR 模型 在 降水 数据 降 尺 度 方面 的 优势 ,基于 MG- 
WR 模型 通过 对 3 个 典型 气候 年 份 CPM 数 据 降 尺 度 
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区 提供 更 高 空间 分 辩 率 降水 数据 ,人 研究 结果 一 定 程 
度 上 可 以 为 我 国 部 分 区 域 相同 降水 条 件 下 的 卫星 
降水 产品 降 尺 度 人 处 理 提供 参考 和 借鉴 ,对 气候 条 件 
多 样 .地 形 复杂 地 区 的 水 文 气象 研究 具有 重要 的 应 
用 价值 。 但 是 ,本 研究 仍 存 在 一 些 问题 有 待 进一步 
研究 。 首 先 ,由 于 研究 区 站 点 数量 较 少量 分 布 不 均 
导致 难以 开展 更 为 细致 的 精度 验证 ,会 降低 结果 的 
准确 度 。 其 次 ,研究 区 域 地 貌 复杂 ,通过 探索 更 多 
其 他 对 降水 变化 敏感 的 环境 变量 ,可 以 建立 一 个 更 
可 靠 、 更 稳定 的 模型 。 


5 结论 


本 文通 过 分 析 MGWR 模型 在 不 同 干 湿 气候 条 
件 下 对 黄河 流域 CPM 卫星 降水 产品 的 降 尺 度 效果 ， 
得 到 主要 结论 如 下 : 

(1) 从 原始 卫星 降水 产品 数据 来 讲 , 黄 河流 域 
2002、2012 年 和 2020 年 GPM 与 地 面 气象 站 点 观测 
值 具有 较 好 的 相关 性 ,GPM 数 据 在 年 际 和 站 点 的 时 
空 尺度 上 精度 较 好 ,在 研究 区 域内 有 一 定 的 适用 性 。 

(2) 利用 MGWR 模型 对 GPM 数据 进行 降 尺 度 
研究 ,其 降 尺 度 结果 与 原始 数据 相 比 ,大 幅度 提高 
了 GPM 数据 在 黄河 流域 地 区 的 空间 分 辩 率 ,能 够 更 
加 准确 地 反映 研究 区 降水 的 空间 分 布 特征 ,并 且 由 
于 MGWR 模型 进行 了 全 局 变量 与 局 部 变量 区 分 ， 
MGWR 数据 较 CGWR 数据 有 着 更 强 的 空间 细节 表现 
能 力 ,在 降水 量 较 丰富 的 2020 年 表现 得 更 加 显著 ， 
说 明 MGWR 模型 把 环境 变量 进行 全 局 变量 与 局 部 
量 区 分 可 以 提升 降 尺 度数 据 的 空间 表现 能 力 。 

(3) 在 黄河 流域 3 个 典型 气候 年 份 中 ,MGWR 
模型 在 2012 年 降水 量 标准 年 中 降 尺 度 结果 表现 较 
好 ,MGWR 数据 比 GWR 数据 的 准确 性 略 有 提高 。 
但 是 MGWR 与 GWR 2 种 数据 精度 在 一 定 程度 上 依 
赖 于 原始 数据 ,其 数据 与 气象 观测 站 点 数据 的 相关 
性 并 未 得 到 提升 ,甚至 损失 了 一 定 的 相关 性 ,这 可 
能 是 由 于 研究 区 域 自然 地 理 条 件 复杂 ,所 选用 的 环 
境 变 量 与 降 尺 度 模 型 无 法 完全 对 降水 的 空间 异 质 
性 进行 解释 造成 的 。 


变 
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Downscaling of GPM satellite precipitation data in the Yellow 
River Basin based on MGWR model 


BAI He", MING Yisen", LIU Qihang", HUANG Chang’ 
(1. Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Xi an 710127, Shaanxi, China; 
2. College of Urban and Environmental Sciences, Northwest University, Xi'an 710127, Shaanxi, China; 3. Institute of Earth 


Surface System and Hazards, College of Urban and Environmental Sciences, Northwest University, Xi'an 710127, Shaanxi, China) 


Abstract: Because the Yellow River Basin of China is a vast area with sparse meteorological stations, limited me- 
teorological data are available. Satellite precipitation data are an alternative for precipitation observations. In this 
study, the precipitation data of the Yellow River Basin for 2002, 2012, and 2020 were considered representative 
of dry, standard, and wet years to downscale global precipitation measurement (GPM) precipitation data. The nor- 
malized difference vegetation index, digital elevation model, slope, land surface temperature, and wind speed that 
reflect the spatial distribution characteristics of precipitation and the characteristics of spatial nonstationarity were 
investigated and used in two downscaling methods, namely the geographically weighted regression model 
(GWR) and mixed geographically weighted regression model (MGWR) to obtain the downscaling precipitation 
data of 1-km spatial resolution in the Yellow River Basin. The downscaling results were verified by the ground 
meteorological station data. The results revealed that: (1) GPM annual precipitation data were highly correlated 
with ground meteorological station observation data in the Yellow River Basin in 2002, 2012, and 2020. (2) Down- 
scaling with the MGWR model considerably improved the spatial resolution. In terms of the spatial details of pre- 
cipitation, the downscaling results of the MGWR model were superior to those of the GWR model. (3) In the 
three typical climate years, the accuracy of MGWR downscaling data in the precipitation standard year was slight- 
ly higher than that of GWR downscaling data. The results of this study can provide a reference for precipitation 
downscaling research in related regions and promote regional climate and hydrological research. 
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